Statistics of Charge Fluctuations in Chaotic Cavities 
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We consider the zero frequency fluctuations of charge inside a mesoscopic conductor in the large 
capacitance limit. In analogy to current counting statistics we derive the characteristic function 
of charge fluctuations in terms of the scattering matrix of the conductor. Using random matrix 
theory we evaluate the characteristic function semi-analytically for chaotic cavities. Our result is 
universal in the sense that it describes not only the fluctuations of charge, but of any observable 
quantity inside the cavity. We discuss equilibrium and non-equilibrium fluctuations and extend our 
theory to the case of contacts with arbitrary transparency. Finally we investigate the suppression 
of fluctuations in the small capacitance limit due to charge screening. 
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The full counting statistics of current fluctuations in 
mesoscopic conductors has attracted the attention of 
many theoretical works during the last decade. In a pi- 
oneering paper Levitov et al. explain the universality 
of current statistics: Coherent charge transfer through a 
two-terminal conductor can be seen as probabilistic pro- 
cess governed by a set of transmission probabilitieai. In 
the following several methods have been developed to 
obtain the full counting statistics of mesoscopic conduc- 
tors: In the original work the authors used a fully quan- 
tum mechanical approach and described the conductor 
by its unitary scattering matrijjSi^. Nazarov provided 
a description in terms of Keldysh-Green's functions that 
can be conveniently applied to conductors with a large 
number of channels^. De Jong characterized double bar- 
rier structures successfully by a fully classical method^ 
based on the exclusion principle. Nagaev proposed a 
semi-classical diagrammatic scheme to obtain higher cu- 
mulants in a systematic mannerii^. Recently, Pilgram et 
al. expressed the full counting statistics in terms of a 
saddle point solution to a semi-classical stochastic path 
integral^. 

In parallel different works addressed the statistics of 
quantities related to current such as the phases in super- 
conducting devicesiSiii or the momentum transfer from 
electrons to a bent conductor—. The interpretation of 
the statistics of these quantities is less intuitive, in gen- 
eral they cannot be characterized by transmission prob- 
abilities only. Momentum transfer for example is not 
quantized in contrast to charge transfer—. 

In this article we report about a quantity of funda- 
mental interest in physics: the fluctuations of charge 
inside a mesoscopic conductor. As soon as parts of 
a mesoscopic system are only coupled capacitively (i.e. 
metallic gates close to a mesoscopic conductor), charge 
fluctuations play an important role: Decoherence for 
instance is described as an average over a fluctuating 
potentialiiii^ which then in turn is coupled electrody- 
namically to the charge. Charge correlations also con- 
tribute to non-equilibrium decoherence when a measure- 
ment process is taking placeiS^. They are as well at 
the origin of Coulomb drag effectsi^. Such models often 
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FIG. 1: Example of a mesoscopic conductor considered in this 
article: the chaotic cavity. We calculate the statistics of the 
charge Q inside the cavity. U is the electrostatic potential of 
the cavity, Vi, V2 are the applied voltages, Ni,N2 the number 
of channels in the contacts, F the contact transparency, kT 
the temperature, and Cg the geometrical capacitance of the 
cavity. 



assume gaussian fluctuations, which is justified in the 
weak coupling limit. When coupling becomes stronger 
non-gaussian fluctuations should be taken into account. 
Methods to calculate higher cumulants are therefore of 
great use! 

In some sense charge fluctuations are diametrically op- 
posed to the fluctuations of current. For a current mea- 
surement we may arbitrarily choose the cross section of 
the conductor, for a charge measurement we are forced to 
define exactly a volume of charge we are interested iniS. 
The choice of this volume will necessarily influence the 
statistics of charge fluctuations. At this point it becomes 
clear that statistics of charge noise are in general non- 
universal in contrast to current noise. There is however 
a class of conductors in which we expect universality to 
be preserved: Disordered systems and in our case chaotic 
cavities. In this case the charge of interest can be char- 
acterized by one single parameter, the dwell time of the 
electrons inside the volume of the cavity2£. Furthermore 
the hamiltonian of a chaotic cavity is random, no basis 
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nor any operator is preferred. Therefore the fluctuations 
of any other internal operator such as the dipolar mo- 
ment will obey the same universal statistics as the charge 
operator. 

In this publication we present a method based on ran- 
dom matrix theorjiSi to calculate the cumulants of charge 
noise in chaotic cavities analytically order by order. For 
the full generating function we show semi-analytical re- 
sults that we confirmed numerically. We emphasize that 
the example of the chaotic cavity in the many chan- 
nel limit can be treated as well in the semi-classical 
framework developed in Ref?2i. As we checked, the semi- 
classical results completely agree with the calculation 
presented here. 

In sectionQlwe derive the generating function of charge 
fluctuations in terms of the scattering matrix. A simi- 
lar expression for the second cumulant has been derived 
long before^S. Results are in particular available for 
the geometry considered in this publication: the chaotic 
cavitjiS^. The charge-charge correlation function has also 
been found at flnite frequenciesSi. These works all em- 
phasize the importance of charge screening, interaction 
effects have to be taken into account! The charge in the 
cavity responds to the fluctuations of the internal elec- 
trostatic potential that is conjugated to the charge. 

We explain in detail the evaluation of the character- 
istic function for an open cavity in equilibrium in sec- 
tion^ We present results for the non-equilibrium case 
and contacts with arbitrary transparency in section ITTll 
We briefly discuss how to include interaction effects into 
our calculation scheme. As an application, we finally es- 
timate the contribution of non-gaussian fluctuations to 
dephasing of electrons passing the chaotic cavity. 

I. GENERATING FUNCTION 

We consider a mesoscopic conductor as shown in Fig. 
n The conductor is described by a hamiltonian Hq . The 
electrostatic potential U on the conductor can be varied. 
The full hamiltonian is therefore H = Hq + UQ where Q 
is the charge operator. In analogy to current statistics^ 
we define the following characteristic function (we set 
;i= 1) 

XtiU) = (e^t(Ho+f Q)e-^t(Ho-f Q)^ 

The first line of Eq. |^ in the Heisenberg picture will 
serve as starting point for our calculation. The sec- 
ond line in the interaction picture contains the Keldysh 
time ordering operator Tk and an integral JJ along the 
Keldysh contour from zero to t and back again. The 
± sign is positive on the upper Keldysh contour and 
negative on the lower contour. The interaction rep- 
resentation allows to verify that the derivative C„ — 
(— lnx(J7)/9?7"|[/=o generates the nth cumulant of 
charge fluctuations. 



The Fourier transform of this characteristic function 
is the probability distribution of charge integrated over 
time /p dt'Q{t'), a quantity without direct physical mean- 
ing. Divided by the elementary charge e it can be under- 
stood as the time spent by all electrons in the cavity after 
time t. Alternatively it can be divided by a geometrical 
capacitance Cg to obtain the probability distribution of 
the phase which is the potential integrated over time t. 
In this publication we mainly concentrate on the charac- 
teristic function itself that generates the moments of zero 
frequency charge fluctuations. 

Using the procedure of Refsi^ the time evolution op- 
erators deflned in Eq. can be expressed by the uni- 
tary scattering matrix of the mesoscopic conductor and 
we obtain 

InXi = 

(2) 

^ / dETi In (1 - HE + nESt(£;, - f )S(£;, +f )) . 

The energy integral reaches from the bottom of the band 
to inflnity. The matrix S(E,U) is the energy and po- 
tential dependent scattering matrix that links incom- 
ing and outgoing current amplitudes. The matrix he 
is diagonal and contains the Fermi occupation factors 
(he)** = n,{E) = [1 + cxp{{E - eVi)/kT)]-^ of channel 
i on its diagonal. In practice it is useful to subtract the 
background charge. We consider separately the charac- 
teristic function at zero temperature and voltage that can 
be expressed by the eigenvalues of the scattering matrix 

^ J° dE i^Y. iMU/2) - 0„(-;7/2))| . 

(3) 

It is obvious from Eq. © that afl even cumulants vanish 
at zero temperature and voltage. The odd cumulants ex- 
cept the first vanish after disorder averaging^^, because 
{d^(l)n/dU^) = for fc > 1. In the presence of disorder, 
subtracting Inxo(t^) from Eq. ^ is thus equivalent to 
subtracting the background charge. We arrive at the fol- 
lowing symmetric expression for the characteristic func- 
tion 

Inx = 

^ dETT In [he + (1 - nE)St (i?, f )S(i?, -f )] 

+ ^ /J" dETiln [(1 - he) + nESt(i?, -f )S(i?, f )] . 

(4) 

At equilibrium one can easily check that Eq. generates 
only even cumulants. The characteristic function then 
only depends on the eigenvalues e*"^" of the scattering 
matrix and is symmetric with respect to U ^ —U . 

II. AVERAGE FOR A CHAOTIC CAVITY 

So far we have been very general. Eq. Q can be ap- 
plied to any mesoscopic conductor, as soon as its scatter- 



ing matrix is knownS^. We now investigate more closely 
the generic example of a many mode chaotic cavity (we 
neglect weak localization corrections). In this case the 
average of products of at most four energy (potential) 
dependent scattering matrices is knownSSiSi. It remains 
nevertheless a non-trivial problem to calculate the av- 
erage of the logarithms in Eq. that contain infinitely 
high powers of scattering matrices. We solve this problem 
in two steps. First we express the characteristic function 
through a Green's function. We then calculate the aver- 
age of this Green's function using a method developed in 
Ref.2i. 

We define the following Green's functions 



r^'?(z,c/) = 



1 



z- A'^'?(t/) 



),F"'=«(z,;7) 



1 



'z - A"'='?(L/)' 



(5) 

The brackets denote the disorder average. The matrices 
A abbreviate products of scattering matrices 



A"'!^ St(-L//2)S(J7/2) 



J^neq^ S 



lr( 



-C//2)Sii(C//2) 



'l2( 



-;7/2)Si2(C//2). 



(6) 

The dimension of the full scattering matrix is iV x , its 
block Sii (the reflection matrix) measures Ni x Ni. We 
introduce the temperature kT, the applied voltage V — 
Vi — V2 > and the integration time t. The characteristic 
functions at equilibrium (kT ^ eV) and for transport 
{kT <C eV) can then be written as integrals over the 
variable z of the Green's functions 



Tr (F'=«(-z, U) + F'=«(-z, -U))) , 



(7) 



and 



^neq(TT\ teV r°o dz 

A V"^ ~" 27r Jo 1 + 2 



(iVi -f (1 -I- z)TrF"'=«(-z, U)) . 



(8) 



To outline the principles of the calculations we first spe- 
cially discuss the simplest case of an open cavity at equi- 
librium. As Brouwer and Biittiker— we introduce the 
potential dependence of the scattering matrix via a vir- 
tual stub described by a M x M reflection matrix 



= Tr*. 



(9) 



The trace (p is linked to the mean level density hy (p — 
27Tdn/{edU) = 2ttNf and to the dwell time by = Gt 
where G is the total dimensionless conductance into the 
cavity and e the elementary charge. The stub is assumed 
to be very large compared to the exits of the cavity; then 
the precise structure of the matrix $ becomes unimpor- 
tant for the calculation. The total scattering matrix of 
the chaotic cavity may be expressed as geometrical sum 

S(C/) = Uaa + Uab (1 - r,(C/)Ubb)"' r,(C/)Uba. (10) 
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FIG. 2: Comparison of semi-analytical and numerical result 
for the logarithm of the characteristic function of equilibrium 
charge fluctuations in an open cavity. We choose M = 160 
for the size of the numerical random hamiltonian describing 
the cavity. 



For an open chaotic cavity the {N + M) x [N + M) ma- 
trix U is distributed according to the circular ensemble 
of unitary matrices. It is divided in blocks Uaa of size 
N X N and Ubb of size M x M. The average over this 
ensemble can be carried out by means of a diagrammatic 
technique^-'^ . The technical details of this average are pre- 
sented in appendix IXI The trace of the Green's function 
turns out to be 



TrF'^ 



iV- 



1 



z-yiz,U) 



(11) 



where y is the root of a cubic Dyson equation 

z~ z{l- ieUr) y-{l + ieUr) y^ + y^ = 0. (12) 

Unfortunately the analytic solution of this polynomial 
equation is rather cumbersome. But we have two ways 
to continue. If we are only interested in the first few 
cumulants of charge fluctuations, we may expand the root 
2/ up to a certain order in the potential U 



y = yo + yiU + y2U' + yy.u^ + 



(13) 



and then solve the cubic equation order by order (As a 
side product of our calculation we can obtain in the same 
way the average Tr((St(-[//2)S(C//2))") by expanding y 
in powers of 1/ z). We then obtain the following expan- 
sion of the characteristic function 

iVi^ {-{eUrf -f ^(e[/r)4 - JL(eC/r)6 + ...]. 

(14) 

The equilibrium fluctuations of the charge in the open 
cavity are non-gaussian. In the next section we will study 



4 



how this result depends on the back-reflection probability 
1 — r of the leads. 

If we are interested in the full characteristic func- 
tion we can find the root of Eq. (|12|l numerically and 
then carry out the integration in Eq. (jT)). The result 
is plotted in Fig. |21 The diamonds indicate a numeri- 
cal check that was performed by numerically averaging 
over a large set of scattering matrices. This fully numer- 
ical solution is very time consuming compared to find- 
ing roots of the Dyson equation. The large U limit of 
the characteristic function can be determined analytically 
to be lnx(?7 — » oo) = — irNtkT / 12. The characteristic 
function % is therefore properly behaved for large times 
t ^ (NkT)^^. As discussed in Refi^ this is the range of 
validity of Eq. (O . Out of equilibrium the applied voltage 
plays the role of the temperature t ^ {NeV)~-^. 

III. RESULTS 

The generalization of our approach to cavities with 
back-reflection at the contacts is not completely straight- 
forward. On the one hand the scattering matrix as given 
in Eq. must be extended and its average (S(C/)) is 
no longer zero. This complication modifies the averag- 
ing procedure for the Green's functions. We discuss the 
crucial points in appendix ^ On the other hand the 
calculation becomes lengthy and requires some algebraic 
software. 

We use the parameters N = Ni + N2 for the total 
number of channels and A = (A^i — A^2)/-^ to describe the 
asymmetry of the cavity (see Fig. The transmission 
probability to pass one of the contacts of the cavity is 
denoted by F. For simplicity we choose it to be the same 
in both contacts. 

For the cumulants of fluctuations at equilibrium we 
obtain 

=iVF^(eT)4(2-r) (15) 

(jeq ^7vr^(eT)6(24- 24r + 8r2). 

The dwell time is given by r and e is the unity charge. 
This result is universal in the sense that it describes also 
the fluctuations of any other operator inside the cav- 
ity. One has simply to replace the constant (er) by 
other units. It is interesting to note that the higher 
cumulants increase with decreasing transparency of the 
leads. We observe the same behavior for non-equilibrium 
charge fluctuations and for current fluctuations in chaotic 
cavities^.. Towards the limit of tunneling contacts the 
non-gaussian fluctuations get stronger. 

In the transport regime the characteristic function can 
be written in the following way 

In^ne, = ^^^jr)L„(A)(zeC/r)". (16) 

n 




FIG. 3: Parameter dependence of non-equilibrium cumulants 
of charge noise. Upper panel: Each cumulant C„ depends 
on transmission F. The prefactors of low cumulants tend to 
zero in the tunneling limit whereas high cumulants stay im- 
portant. Lower panel: Each cumulant depends on the asym- 
metry X — {Ni — N2)/{Ni + N2). Cumulants higher than two 
may change their sign. 



The functions Kn and L„ describe the dependence of the 
cumulants on transparency of the contacts and asymme- 
try of the cavity. The first few coefficients are given by 

Ki=T if2 = r(2-F) ifa = F(3 - 3F + F2) 

r _ l+A T _ l-A^ T _ A(l-A") 
i^l - — ^2 - — g— ^3 - 12 

(17) 

These factors are plotted in Fig. as a function of trans- 
parency and asymmetry. The higher the cumulant, the 
slower the factors Kn diminish towards the tunneling 
limit. The plot indicates a non- analytic behavior for 
F and n — s- 00. A similar effect occurs in extremely 
asymmetric cavities at |A| ~ 1: Non-gaussian fiuctua- 
tions are enhanced. 

Up to now we considered a completely non-interacting 
system {Cg ^ e'^Np). We now relax the condition 
Cg ^ e^Np by including the effect of charge screen- 
ing. In the many-channel limit considered here Coulomb 
blockade effects are unimportant. Interaction effects can 
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then be treated by including a self-consistent response of 
the electrostatic potential U (the bottom of the band) 
to the charge fluctuations^'^. We introduce a fluctuating 
electrostatic potential U — Uo + SU (t) that is coupled ca- 
pacitively to the fluctuations of the total screened charge 
on the cavity 6Qscr = Cg6U. We then have to include the 
response of charge cuniulants C„ to fluctuations of the 
potential 5U, the so called cascade correctiongi^. In the 
case of cavities this procedure is almost trivial, because 
only the first cumulant of charge Ci, the mean charge, 
depends on the electrostatic potential, see appendix [HI 
It turns out that the screened cumulants differ from the 
unscreened cumulants by a universal prcfactor (in the 
linear bias regime) 



1 



(18) 



An identical prefactor has been obtained in Ref.-- for 
the second cumulant. This simple result relies on two 
facts: All cumulants higher than one do not depend on 
the electrostatic potential, because correlators of scatter- 
ing matrices {{S'^ {Ui)S{U2))'^) only depend on the differ- 
ence U1 — U2 after disorder averaging. Furthermore, the 
density of states Np depends in principle on the parti- 
cular scattering matrix. Only in the many channel limit 
we can consider Np and C„ as being independent when 
averaging over disorder. In general a relation as simple 
as Eq. ifTFI) cannot be expected! 



IV. APPLICATION: PHASEFLUCTUATIONS 

In this chapter we estimate the contribution of non- 
gaussian fluctuations to the dephasing of electrons pass- 
ing a chaotic cavity. This allows us to formulate a validi- 
ty condition for the gaussian approximation. Models for 
dephasing in chaotic cavities have been considered for 
instance in Refsi^SiSS. Dephasing rates in cavities have 
been measured by Huibers et al.™ and Hackens et st\M>. 
Dephasing due to charge fluctuations has been treated in 
the scattering formalism in Refsi^S'SSiS. 

We approximate the dynamical phase picked up by an 
electron in the cavity by 



~ e 



r dt'u{t') r dt'Q'-'-it'). (19) 

Jo Jo 



The screened charge is denoted by Q'"^'', its statistics is 
in the long time limit given by Eq. H18(l . We consider 
the phase (/) to be a classically fluctuating fieldi^ (with 
non-gaussian fluctuations) and average the phase over 
equilibrium charge fluctuations 



scr I __ r^scr 



(20) 



Here we assumed additionally that fluctuations on the 
time scale of r behave as fluctuations in the long time 



limit. This is justified at high temperatures kT 
if the frequency dependence of the response function 
d{Q)/dU (see appendix Q may be neglected. We com- 
bine Eqs. H15|) and l|18|) . and insert the cumulants into Eq. 
(|20|1 . A comparison of second and forth cumulant gives 
the following condition for the validity of the gaussian 
approximation 



Cl2 



r 



G2 C2 



12 



< I 



(21) 



where we introduced the electrochemical capacitance 
= Cg"^ + [e^Np)'^ and the dimensionless conduc- 
tance G — TN . Note that < Cg, therefore condition 
(|21|l always holds in the large conductance limit consid- 
ered in this article. Nevertheless, Eq. I|21|) indicates that 
non-gaussian corrections become important in the limit 
of few channels. They are strongest in the charge neutral 
limit Cn = Cg and increase in the tunneling limit F — > 0. 



V. CONCLUSIONS 

In this article we described a way to calculate higher 
cumulants of charge fluctuations inside a mesoscopic con- 
ductor in the long-time limit. We first derived a general 
generating functional for such fluctuations in terms of the 
scattering matrix of a non-interacting mesoscopic system. 
As a specific example we considered a chaotic cavity and 
showed how to average the generating functional over dis- 
order using random matrix theory. We then studied the 
dependence of the generating functional on parameters 
as the asymmetry or transparency of the leads attached 
to the cavity. We find that higher cumulants contribute 
most to fluctuations in the tunneling limit. Finally we 
calculated the suppression of fluctuations in the pres- 
ence of screening. Although higher cumulants of charge 
fluctuations are probably not directly measurable, their 
knowledge is of great practical importance. It allows for 
instance to study decoherence effects due to fluctuations 
beyond the gaussian theory. 

The authors thank P. Samuelsson and E. V. Sukho- 
rukov for valuable discussions and acknowledge the sup- 
port of the Swiss National Science Foundation. 



APPENDIX A: DISORDER AVERAGE OF 
GREEN'S FUNCTIONS 

In this appendix we discuss how to average the po- 
tential dependent Green's functions jSJl in the equilib- 
rium case. In Ref.^^ Brouwer and Becnakker present a 
procedure to calculate the distribution of transmission 
eigenvalues of chaotic cavities from a disorder averaged 
Green's function. We adapt this procedure to our pur- 
poses. In a first step the series in Eq. (|10|) is reexpressed 
in a compact formulation by introducing a set of matrices 
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S(C//2) 
SH-U/2) 



F(z) = 



F'(z) 
F(z) 



T\ R' 



r' 






T' 








1 

1 



t 
tt 






(r')t 




rH-U/2)/ 



(Al) 

The transmission and reflection matrices t, r and r' de- 
scribe the back-reflection in the contacts. Without bar- 
riers we have t = 1 and r = r' = 0. The fluctuating part 
of the scattering matrix then is 



(5S = T' ( 1 - UR' 



U 



U 

ut 



(A2) 



where U is the uniformly distributed unitary random ma- 
trix. The matrix Green's function F that contains the 
desired Green's function (O as off-diagonal element can 
be written as difference F = F+ — F^. Both components 



1 



1 



2^ ViipSC 



(A3) 



can be averaged separately. According to Brouwer and 
Beenakker^^ they obey a Dyson equation 



(A4) 



F± =X± 1 + S±(F±)F^ 



with a self-energy matrix = S^(F^) depending on 
the Green's function F^. The matrix is defined as 
X = R' 4- TCT'z-i/2. The equation for the self-energy 
may now be expanded up to first order in 1/M where 
M is the dimension of the virtual potential dependent 
stub 0. After impurity averaging and some algebra we 
obtain that the self-energy is of the form 



1 



(A5) 



and fulfills Eq. (P|l . 



In order to obtain the density of transmission eigenvalues 
as in Refi^ one chooses 



Ci 













2 







(B2) 



and exploits the fact that Ci,C2,S commute and 
C1C2 = 0. These conditions do not hold in our case. 
The matrix C contains more entries. In equilibrium for 
instance we need Ci = C2 = 1- We overcome this prob- 
lem by replacing C through a matrix C that encounters 
the back-refiection in the leads 



c 



C 1 RCz 



-1/2 



R 







(B3) 



One can easily verify that this substitution generalizes 
the Green's function introduced in Eq. (|A3|I to the case 
of barriers at the contacts. The rest of the calculation 
follows the lines of appendix 1X1 



APPENDIX C: SCREENED CHARGE 
FLUCTUATIONS 



In this appendix we explain the derivation of Eq. I|18l) 
using the cascade principle (for an introduction to cas- 
cade corrections we refer the reader to Refsii'2*^). Charge 
and potential fluctuations of the cavity are linked via the 
capacitance Cg: 



SQs 



CgSU, SQs 



(CI) 



The charge fluctuations are composed of bare fluctua- 
tions at constant potential SQ and the linear response 
d{Q)/dU — —e^Np to 5U (the screening charge). The 
fluctuations of 6Q are known from the non-interacting 
problem (see Eqs. ((T3|) and (fTE|) '). Solving Eq. HC1() for 
the total charge fluctuations SQscr the third cumulant 
without cascade corrections may be written as 



APPENDIX B: EXTENSION TO CAVITIES 
WITH BARRIERS 



{{SQscrf) = 1 + 



~c7 



(C2) 



In the presence of back-reflection at the contacts we 
have to modify the summation carried out by Brouwer 
and Beenakker^^. Now the scattering matrix of the leads 
is given by t = ivT and r = r' = \J\ — V. The additional 
complication due to barriers arises from the fact that the 
mean of the scattering matrix S is no longer zero. The 
needed form of the matrix C is 



The two possible cascade correction are given bjii 

3 {bUdQscr), 3 — — [dUdQscr) ■ 

(C3) 

But these corrections are zero, because {{SQscrY) does 
not depend on U and (SQscr) depends linearly on U. 
Cascade corrections to all higher cumulants vanish in the 
same way. Eq. (|18|) is therefore valid for arbitrary n. 
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